Electronic properties of the M0S2— WS2 heterojunction 
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We study the electronic structure of a heterojunction made of two monolayers of M0S2 and 
WS2. Our first-principles density functional calculations show that, unlike in the homogeneous 
bilayers, the heterojunction has an optically active band-gap, smaller than the ones of M0S2 and WS2 
single layers. We find that that the optically active states of the maximum valence and minimum 
conduction bands are localized on opposite monolayers, and thus the lowest energy electron-holes 
pairs are spatially separated. Our findings portrait the M0S2— WS2 bilayer as a prototypical example 
for band-gap engineering of atomically thin two-dimensional semiconducting heterostructures. 



Engineering the electronic properties of semiconduc- 
tors by using hetero junctions has been the central con- 
cept in semiconductor science and technology for 5 
decades With the advent of quantum wells, band- 
gap engineering of quasi- two dimensional semiconductors 
made it possible to observe a wealth of new physical phe- 
nomena, including the integer and fractional quantum 
Hall effects in modulation doping GaAs/GaAlAs^i^ the 
condensation of both excitons in double GaAs quantum 
wells of GaAs,^ and exciton-polaritons in II- VI quantum 
wells^ii and, more recently, the Quantum Spin Hall phase 
in CdTe/HgTe quantum wells. ^ 

The isolation^ of truly two dimensional crystals, such 
as graphene and M0S2, and their use to fabricate field ef- 
fect transistors ^iii has opened a wealth of new venues in 
physics and material science in general, and more specif- 
ically in the design heterostructures Thus, graphene 
bilayer^ and graphene on boron nitride^ have different 
electronic properties than freestanding graphene. 

The properties of bulk M0S2 and its nanostructures, 
such as nanotubes^ii^ fullerenes^ and nanoislands,^^ 
have been studied for a long while, including even chem- 
ically exfoliated single planes. More recently, the study 
of electronic and optoelectronic devices based on a single 
M0S2 layer has taken impetus for several reasons. First, 
it was found that monolayers of M0S2 feature a direct 
band gap of 1.8 eV with strong photoluminescence^i^ 
as opposed to bulk M0S2 which has indirect band gap 
of 1.29 eV. Second, the fabrication of a high mobility 
field effect transistor based on single M0S2 layer has 
been reported^^^ Third, the combination of hexagonal 
symmetry, large spin-orbit coupling and lack of inversion 
symmetry, give rise to gapped graphene like bands, with 
two valleys and strong spin- valley coupling Taking ad- 
vantage of these unique properties, optical spin pump- 
ing is turned into valley-polarized photo carriers^ — 
which opens new possibilities in the emerging field of 
valley tronics.^^ 

Importantly, other transition metal dichalcoghenides, 
such as WS2, as well as MoSe2 and WSe2 are expected 
to have similar properties^ — and the first experimen- 
tal demonstrations of monolayer WS2 have just been 
reported.^- All of the above naturally leads us to investi- 
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FIG. 1. (Color online) Schematic views of the M0S2-WS2 
heterojunction of different stacking (i.e. 07, 027, AA, T). 
Each stacking is obtained either by a monolayers translation 
T and /or a rotation R with respect to each other. Red, gray, 
and yellow spheres represent W, Mo, and S atoms respectively. 



gate the electronic properties of transition metal dichal- 
chogenide (TMD) multilayer^—. Here we report our 
results on the simplest case, a bilayer of M0S2 and WS2, 
which both have the same crystal structure and very sim- 
ilar lattice constant. In particular we are interested in 
how the stacking of different TMD monolayers (see Fig. 
[1]) can result in heterostructures with electronic prop- 
erties different from the homogeneous TMD monolayers 
and multilayers. 

Our calculations were performed with the Vienna ab- 
initio package (VASP)^ based on the local density- 
functional approximation^ plane-wave basis (£^cut = 
400 eV) and non-collinear projector-augmented waves 
(PAW) method^S;^!^ We treat the both the transition 
metal orbitals 4p, 55, Ad together with the Sulphur orb- 
tials 3s and 3p as valence states, and the rest are con- 
sidered as core. We use the Perdew-Burke-Ernzerhof's^- 
version of generalized gradient approximation to describe 
the exchange correlation density functional. All calcu- 
lations are carried out using a 1x1 supercell with vac- 
uum thickness not smaller than 17 A. The F-centered 
Monkhorst- Pack's^ mesh (6 x 6 x 1) of the k-points was 
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FIG. 2. (Color online) Band structures of: a) M0S2 mono- 
layer, b) WS2 monolayer, d) M0S2 bilayer, e) WS2 bilayer, 
f) M0S2— WS2 heteroj unction. The stacking of bilayers is C7 
(see Fig. [Hb)). c) scheme of the BZ with the line along which 
the band structures are calculated. Evac stands for vacuum 
energy. The Fermi energy lies at the intersection of white and 
yellow regions. 



used to sample the BZ. 

For reference, we discuss first the electronic properties 
of isolated M0S2 and WS2 monolayers (MLs)^-'^^^ 
The crystal structure of 2H-M0S2 (2H-WS2) consists of 
two 2D parallel triangular lattices of S atoms separated 
by the same lattice of Mo ( W) atoms translated by 1 /3 
of the unit-cell diagonal, with lattice constant a = 3.19 A 
(a = 3.20A).^^ The corresponding Brillouin zone (BZ) is 
also hexagonal, with two inequivalent K and K' points 
(valleys). We show the corresponding energy bands in 
Figs. [2fa) and[2](b)), which are in agreement with pre- 
vious work using the same methodology j^^i^^ Both MLs 
are direct band semiconductors with a maximum valence 
band (VB) and minimum conduction band (CB) located 
in the K and valleys. The band-gap values equal 
1.58 and 1.50 eV for M0S2 and WS2 respectively. Our 
calculations also show that, when referred with respect 
to the vacuum energy, the band structures of both MLs 
are shifted (cf. Figs. [2](a) with[2](b)), on^ account of the 
different electronegativity of the Mo and W. 

Because of the lack of inversion symmetry and a strong 
spin-orbit coupling (SOC) the valence and conduction 
bands are spin-split at the K and points. The sign 
of the spin splitting changes from K to resulting in 
the so called strong spin-valley coupling;^ The splitting 
is higher in WS2 ML (435 and 27 meV for VB and CB 
respectively) than M0S2 ML (147 and 3 meV for VB and 
CB respectively) due to the higher atomic number of W 
than Mo. 

The band-dependence of the spin splittings is ac- 
counted for by the atomic orbital composition of the 
states. Our population analysis reveals that, at the 
K point, the CB minimum is mostly made by Mo 



(/ = 2,m = 0) orbitals, whereas the VB maximum 
dominant contribution comes from the dxy and dx'2-y2 
{I = 2,m = ±2) orbitals. To leading order in the 
SOC, this should yield a non-zero valley dependent spin- 
splitting only in the VB, in agreement with the toy model 
proposed by Xiao et al^ However, both VB and CB 
states at the i^, K' points have small contributions com- 
ing from the Sulphur px and Py orbitals (/ = 1, m = ±1). 
These are probably behind the small but finite splitting 
in the CB. 

We now discuss the electronic properties of the bilayers 
that can be formed stacking the WS2 and M0S2 mono- 
layers. We have verified that the main features of the 
electronic structure are quite insensitive to the stacking 
(see Fig. [T]for the different stackings), thereby we focus 
on the band structure of the C7 stacking (see Fig. [Ifb)) 
presented in Fig. EJf). This is the stacking of bulk M0S2 
and WS2. Comparison of monolayer and bilayer bands 
in Fig. 2 indicates that interlayer coupling is not strong. 

The electronic structure of Mo-Mo and W-W bilayers 
(Figs. EKd), (e)) can be rationalized in terms of two con- 
cepts: interlayer coupling of degenerate monolayer states, 
which splits most of the monolayer states, and the exis- 
tence of a symmetry center in the C7 stacking, which 
prevents spin splittings. The interlayer splitting is signif- 
icantly stronger for the valence band at the F point than 
for the VB and CB at K points. As a result, the highest 
VB state moves to the F point for the W-W and Mo-Mo 
bilayers, which become indirect gap systems. 

In the case of the Mo-W heteroj unction the inter- 
layer coupling competes with the energy difference of the 
monolayer states, shown in Fig. 2 (a),(b). As a result, 
the VB at the F point is almost degenerate (A£^vb = 27 
meV) with the top of the VB at the K and K' points. 
Consequently, a significant population of photoexcited 
holes will be available at the K and K' points, and pho- 
toluminescence will be not quenched. In this sense, the 
M0S2— WS2 heteroj unction - unlike the homogeneous bi- 
layers - will be optically active. In addition, the Mo-W 
bilayer does not have inversion symmetry, so that spin 
splittings at the K and K' points occur like in the mono- 
layers. 

A summary of the electronic states in the neighbor- 
hood of the K point, relevant for the inter-band optical 
experiments, both for M0S2, WS2 and the Mo-W bilayer 
is shown in Fig. [3l It is apparent that interlayer cou- 
pling at this point is negligible and the bilayer bands are 
nothing but a superposition of the monolayer states. As 
a result, the top of the VB is in the W layer and the 
bottom of the CB name lies in the Mo layer, forming 
a type II structure.^ In addition the Mo-W bilayer gap 
is 1.2eV, 0.3eV smaller than the gap of the monolayers. 
In contrast, the top of the VB at the F point is delo- 
calized in both planes. The resulting scheme of levels is 
shown in Fig. 3(d). We expect that intra-layer transi- 
tions have a stronger quantum yield, on account of their 
larger electron- hole overlap,^ but relaxation to the lower 
energy spatially separated electron-hole pair is expected. 
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FIG. 3. (Color online) Zoom of band structures at the K point for : a) M0S2 monolayer, b) WS2 monolayer, c) M0S2— WS2 
bilayer, d) scheme of possible optical transitions in the M0S2— WS2 bilayer. Blue (red) lines describe the bands of the states 
localized on the W (Mo) atoms. Bold (dashed) lines describe the states of spin up (down). 



We finally comment on the limitations of our method- 
ology. First, PEE is known to underestimate the band 
gap. The use of either non-local functionals^^ and/or 
GW approximation^!^ yields a better agreement with 
the experiments. Second, excitonic effects, not included 
in the band structure calculation, are known to produce a 
large shift in the absorption threshold and photolumines- 
cence peaks. In spite of this, we expect that the band- 
gap reduction and the segregation of electrons and holes 
in different atomic planes will be confirmed by experi- 
mental work and/or more sophisticated methodologies. 

In summary, we have studied the electronic proper- 
ties of the M0S2— WS2 system as an example of transi- 
tion metal dichalcogenide two-dimensional heterostruc- 
ture. We find that, in contrast to the Mo-Mo and W-W 



bilayer s, the band-gap is direct. Additionally we find 
that the lowest energy electron and highest energy hole 
states in the optically active K point are localized on dif- 
ferent monolayers. In this sense, the Mo-W bilayer forms 
a type II heterostructure. The combination of band gap 
engineering in hetero junctions found here, together with 
the reported electrical control of electronic and optical 
properties in these system o^^i^^i^^i^^ , hold the promise 
of a bright future for optospintronics in two dimensional 
transition metal dichalcogenides. 
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